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Abstract. A new class of analytical 2-D solutions of the full set of the steady magnetohydrodynamic (MHD) 
equations, describing an axisymmetric helicoidal magnetized outflow originating from a rotating central object, is 
presented. The solutions are systematically obtained via a nonlinear separation of the variables in the momentum 
equation. The analysis yields three parameters which measure the anisotropy in the latitudinal distribution of 
various flow quantities. Topologically, the wind speed is controlled by an X-type critical point that acts to filter 
out a single wind-type branch and the Alfven singularity. The solutions can be regarded as an extension outside 
the equatorial plane of the Weber & Davis (1967) model of magnetized winds but with a variable polytropic index. 



Key words. Magnetohydrodynamics (MHD) 
- ISM: jets and outflows 



plasmas - Sun: magnetic fields - solar wind - stars: winds, outfiows 



1. Introduction 

Recent detailed observations of the solar wind (e.g., via 
Ulysses and SoHO, Feldman et al. 1996) and other astro- 
physical outflows (e.g., via the Hubble Space Telescope, 
Biretta 1996) have highlighted the need to provide a satis- 
factory description of cosmic MHD outflows which deals 
adequately and simultaneously with their: (i) spatial 2-D 
(or 3-D) character, (ii) nonlinear dynamics, (iii) detailed 
energy balance and (iv) time-dependent nature. At the 
moment this task is far from being completed and still re- 
mains a challenge for the future. The present treatment is 
not going to change this fact that we do not have available 
to this day a satisfactory MHD solar/stellar wind descrip- 
tion. Hence, by necessity, in the present paper we shall 
constrain our efforts only to the above items (i)-(ii). 

The first models of astrophysical outflows were poly- 
tropic and 1-D (Parker 1958, 1963), or, all forces per- 
pendicular to the equatorial plane were ignored (Weber 
& Davis 1967). The discovery of flows through coronal 
holes called for the investigation of quasi-two dimensional 
polytropic models at flrst (Kopp & Holzer 1976, Habbal 
& Tsinganos 1983) and fully 2-D numerical models later 
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(Pneuman & Kopp 1971, Steinolfson et al. 1982, Nerney 
& Suess 1975, Sakurai 1985). Recently, time-dependent 
studies of thermally and/or magneto-centrifugally-driven 
polytropic wind-type outflows provided new and promi- 
sing results in relation to the question of the degree of col- 
limation of the outflow (Ouyed & Pudritz 1997, Keppens 
& Goedbloed 1999, Ustyugova et al. 1999, Tsinganos & 
Bogovalov 2000, Krasnopolsky et al. 2000). Non-polytropic 
wind modelling on the other hand, with energy and mo- 
mentum addition and finite thermal conductivity has been 
increasingly used because observations have highlighted 
the fact that the acceleration of the solar wind in high- 
speed streams (Feldman et al. 1996, Giordano et al. 2000) 
does imply energy and momentum addition in the solar 
corona (Leer & Holzcr 1980, Steinolfson 1988, Suess et al, 
1996, Hansteen et al. 1997, Wang et al. 1998). 

A fully 2-D MHD modelling of steady hydromagnetic 
wind-type outfiows in open magnetic fields has been in- 
troduced (Low & Tsinganos 1986, Tsinganos & Low 1989) 
with a variable polytropic index 7 wherein it was shown 
that the density needs to decrease with latitude for an ac- 
celerated wind. This is because a dipolar magnetic field 
needs to be kept open by a pressure that must decrease 
towards the pole. If the density does not vary with lati- 
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tude, there is a smaller pressure gradient to drive the flow 
near the pole, exactly where the magnetic field is open to 
allow the wind to escape. The resulting acceleration is too 
low, since gravity dominates, and the flow does not reach a 
high enough terminal speed. The only way out is to allow 
the density to increase with latitude, faster than the pres- 
sure does (Hu & Low 1989). Subsequently, solutions of 
the MHD equations with a latitudinally-dependent den- 
sity and a helicoidal geometry of the streamlines were 
analysed (Tsinganos & Trussoni 1991, hereafter referred 
to as TT91), but the latitudinal dependence of the dif- 
ferent quantities was assumed a priori. Under the same 
assumptions, this work was later continued (Trussoni & 
Tsinganos 1993, Trussoni et al. 1997). Other analytical 
solutions which have used similar assumptions have con- 
centrated on the asymptotic analysis and its links to the 
intitial boundary conditions and were able to derive a 
criterion for the coUimation of winds into jets (Sauty & 
Tsinganos 1994, Sauty et al. 1999). 

In this paper we shall deduce from the governing MHD 
equations via a nonlinear separation of the variables the 
latitudinal variation of the density and other relevant 
physical quantities of the wind, instead of adopting a pri- 
ori for them a specific form, as it was done in TT91 and 
in subsequent papers of that series. This shall also fur- 
ther extend previous work (Lima & Priest 1993, hereafter 
referred as Paper I), which deduced a general class of so- 
lutions of the hydrodynamic equations relevant to stellar 
winds with a helicoidal geometry. Sect. 2 is devoted to the 
method of solution. Sect. 3 deals with its parametric study 
and finally Sect. 4 ends with a discussion of the results. 

2. Method of solution 

2.1. Basic equations and assumptions 

We shall seek 2-D solutions of the set of steady MHD 
equations describing the dynamical interaction of an in- 
viscid, compressible and highly conducting plasma with 
an axially symmetric magnetic field created by a central 
rotating object. These include Maxwell's equation for the 
divergence of the magnetic field, the induction equation 
in the limit of large conductivity and the equations for 
conservation of mass and momentum (Priest 1982). The 
pressure is related to the density and temperature through 
the classical ideal gas law p = {2kB/m)pT while the tem- 
perature is obtained from an energy conservation equation 
such as the first law of thermodynamics, including a heat- 
ing/cooling term. 

The system of MHD equations can be closed if the 
heating function is functionally related to the thermody- 
namic variables and the fiow field. One such simple exam- 
ple which has been widely used in the literature is when 
we have a polytropic law between pressure and density. 
The alternative approach that we follow here is to find 
some other functional relationship between the heating 
rate, pressure, density and fiow speed such that the varia- 
bles are separable. 



The obvious choice of coordinate system appropriate to 
this type of problem is that of spherical polar coordinates 
(r, 0, (/)), with 9 as the co-latitude. With axisymmetry, in 
order to make the above equations more tractable from 
an analytical point of view, the simplifying assumption 
will be made that the meridional components of both the 
velocity and magnetic field can be neglected (i.e. Ve = 
Bg = 0). Under this assumption the projection of the 
field lines on the poloidal plane corresponds to radial lines. 
Note that it has been recently shown that this assumption 
of radiality is a good one, at least in the case of the solar 
wind (Wang et al 1998, Tsinganos & Bogovalov 2000). 

Then, the system of MHD equations can be re- written 
explicitly in the following way: 



d 



(By) = 0, 



— {rV^Br-rVrB^)^0, 
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dr 
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2.2. A technique based on a separation of the variables 

At this stage, we introduce a second simplifying assump- 
tion, namely that the variables are separable in r and 9. 
This will transform the above system of partial differential 
equations into a system of ordinary differential equations 
which are analytically more tractable. 

Denoting by rg the radius at the base of the atmo- 
sphere, we can non-dimensionalise all quantities with re- 
spect to their values at this reference level. In particular 
we set R = r/rp. Using the assumption of separation of 
variables, we can write the radial velocity, the azimuthal 
velocity and the azimuthal magnetic field as, respectively. 



Vr{R,9)=VoY{R)vr{9), 
V4R,9) = ViA{R)v49), 
B^{R,9) ^ BiMiR)b^{9), 



(7) 
(8) 

(9) 



where Vq, Vi and Bi correspond to their reference values. 
Note that, at this stage, the functions Vr{9), v^{9) and 
&0(0) are completely arbitrary. /^From Eqs.(|^), (^) we must 
also have 



Br{R, 9) — Bq 



br{9) 
i?2 ■ 



p{R,9)^Po 



(10) 



with Bq and po as the values of the radial magnetic field 
and density at the reference level, while hr{9) and g{9) are 
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as yet to be determined. Note that the radial magnetic 
field has a monopole geometry modified by the presence 

of brie). 

Using the same technique described in Paper I, let us 
eliminate the pressure term between the r- and 9- com- 
ponents of the momentum equation, Eqs.lQ), (|^), by dif- 
ferentiating the first one with respect to 6 and the second 
one with respect to r, and then adding. The resulting ex- 
pression is 

89 r d9V'"dr in 89 \ r 



'An 09 ' ' 



d ( 2 cos 



sitl9 



dB£\ _ d_ f pGM 
'~d^J ^ 09 \ r2 



4TTdr\ ''09 



(11) 



J_0_ cos 9 
An Or \ sin 9 



J__9 

An Or 



B, 



OB^ 
' 09 



Under the assumption of separation of variables, the above 
equation will be transformed into an ordinary differential 
equation involving functions of R alone. For that purpose, 
the functions of 9 in each term will be set proportional to 
each other. We should keep in mind, however, that this 
last procedure may not give us the most general separable 
solution in the sense that some terms could be also set 
proportional to the sum of the others. Thus, from the last 
two terms, we may write 



dbd. 



ebi 



, cos 9 



^ d9 sm f 

which implies that 

■ bl ^ 



_d 
d^ 



sm 



2e , 



= 0, 



and so 



sm 



(12) 



(13) 



(14) 



where e is an arbitrary constant. Note that the constant of 
integration has been set to unity, without loss of general- 
ity. Any other value for this constant can be incorporated 
into the radial dependence of B^{R,9). In what follows, 
the choice of the constants of proportionality between dif- 
ferent terms as well as the constants of integration will be 
made so as to obtain the simplest possible solutions, with- 
out loss of generality. Comparing the 7*'* and 8*'* terms, 
we can put 



dbr 



br 

which implies 



, cos 9 



9 

d^ - ^'^^^^ 



br{9) = \Jl + Msm^'fi*, 

where is a second arbitrary constant. 



(15) 



(16) 



From the 5*'* and 
obtain 

dg 2COS0 
a9 ^ sm 9 

thus resulting in 

5(61) = l-f(5sin2<^6', 



3th 



terms in Eq.(pi]) we may also 

(17) 

(18) 



in which 5 is the third arbitrary constant. As for the 6 
and 8*'* terms we can write 

9^1 = hi, 
or, equivalently, 

sin' 6* 



th 



(19) 



MO) 



VT+Tsin^ 



(20) 



Finally, from the (/)-component of the momentum equa- 
tion, Eq.(^), for the separation in R and 9 to work we 
must have 



gVrV^ = brb^, 



givmg 



Vri9) 



'l + nsm^^9 
l + Ssm'^'9' 



(21) 
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Fig. 1. Co-latitudinal dependence of the radial magnetic 
field br{d) (solid line) and of the radial velocity (dashed 
line) for e = 1 and 6 ~ A. 

Note that from the four remaining terms in Eq. ( pT| ) — 
I'**, 2"^^ 3'''^ and 4*'' — we obtain identities, provided the 
proportionality constants are chosen accordingly. 

We have thus deduced general ^-dependences for the 
density and hydromagnetic field, under the assumption of 
separation of variables, Eqs.(|l|), (|l|), ([8|), (|o|) and (H), 
in terms of three arbitrary constants: {d,e,fj,). These pa- 
rameters effectively control the anisotropy in the outflow, 
in the magnetic field and in the density distribution. 

First, the parameter S is related to the ratio of the 
density at the equator {9 = n/2) to that at the pole {9 = 
0) and so the higher it is, the more the density distribution 
deviates from the spherically symmetric case (S = 0). 
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Second, the parameter fi is related to the ratio of the 
radial kinetic energy density at the equator to that at the 
pole. Finally, the parameter e controls the width of the 
profile of the speed and density for some fixed variation 
between pole and equator. 

Note that in the TT91 model, the simplest possible 
forms of br{0), g{9) and Vr{9) were chosen a priori such 
that they were able to simulate existing observations of 
the solar wind. Their expressions constitute a special case 
of the forms deduced in this work, corresponding to e = 1 
and fj, = —1. 



(a) 



2.3. Angular momentum and angular velocity 

Returning to the 0-component of the momentum equa- 
tion, Eq.(^), it can be manipulated to give 



d_ 

dr 



r sin OVj, — r sin ( 



^ BrB^ 



0. 



(23) 



The solution of this equation introduces a function of of 
the form 



L{e) = rs\TLeV^-rs\TL0- 



(24) 



which is the total angular momentum per unit mass 
loss carried away by the wind, along each flow line 
9 = const (see Weber & Davis 1967). The first term 



on the right-hand side of Eq.(24) is the angular momen- 
tum carried by the advection of the flow while the second 
term represents the torque associated with the magnetic 
stresses. Analogously, the solution to the induction equa- 
tion, Eq.(^), introduces another function of 9 of the form 



m = ^Av^-B^^ 

r sm 9 \ Br 



(25) 



which corresponds to the angular velocity of the roots of 
the field lines on the surface of the central object. 

The azimuthal components of the velocity and magne- 
tic fields can now be obtained via these two functions L{9) 
and ri{9). After some straightforward manipulation we ar- 
rive at: 



V - 



X^-J^]Rsin9m, 



Bs — 



Boro (R^/Rl~l\ n{9) sin 9 \/ 1 + S sin^" 9 



VoY, 



1 



Ml 



R 



(26) 



(27) 



where = {Vt/Va)'^ is the radial Alfven Mach number 
which corresponds to the ratio of the radial velocity, Vr, to 
the radial Alfvenic velocity, Va — Br/y/inp, and (i?*, F*) 
is the point for which AIa = 1. At this point, Y^,R^ = 



wher 



^0 



(28) 



is the ratio of the radial velocity to the Alfvenic radial 
velocity at the base of the wind. A physical interpre- 
tation of the radial Alfven Mach number Ma can be 




Fig. 2. Co-latitudinal dependence of the angular velocity 
of the roots of the field lines il.{9): in (a) for 5 = 4, in (b) 
for e — 2, in (c) for e = 10 



found if we write Ma = {pV^ /2)/{Bf/8'K), which rep- 
resents the ratio of kinetic to magnetic energy density. If 
Ma <C 1 magnetic energy dominates, whereas if Ma ^ 1 
the dominant energy is kinetic. The radius i?, represents 
the distance at which magnetic energy ceases to dominate 
over the kinetic energy. The regularity condition at i?*, 
Ma{R*) — 1, constrains the two functions L{9) and ^l{9): 
L{9) ^ n[9)rl SIT? 9. 

If we now substitute Eqs.(|2^), (^ into Eqs.(|), (|), 
respectively, and compare the resulting expressions with 
Eqs.(|6|) and (|2^), we can deduce that VqBi = ViBo, to- 
gether with the general form for ^{9), 



n{9) 



AVqK, sin 



'^0 y/l+J^n^ 



(29) 
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where A — Vi/Vq = Bi/Bo measures the ratio of az- 
imuthal to radial velocities at the base of the wind, while 

m ^ (30, 



+ a sinks' 

The angular dependence of n{9) as the parameters e and 
S vary is worth some attention and is plotted in Fig. ^. 
This behaviour of ^{0; e, S) which emerges naturally from 
the separation of the variables in the governing equations 
can be compared, for example, with the rotation law of 
sunspots and solar photospheric magnetic fields. Except 
for the regions around the poles (since from Eq. (P9|), 
n{9) goes to zero for 9 = and 180°), the graphs of Figs. 
^(a,b,c) for e > 1 show similarities with the angular de- 
pendence of r2(0) arising from the yearly averaged rotation 
profiles of photospheric magnetic fields (Snodgrass 1983). 

2.4. Mass and angular momentum efflux 

A quantity of interest is the mass efflux (pV^r^) or, equiv- 
alently, the mass loss rate per infinitesimal solid angle dE 
at the angle 9 (see TT91). From the equation of conser- 
vation of mass Eq.(^) this has to be a function of 9 alone, 
which we shall denote by m{9) 



m{0) = PoVorl^{l + fi sin^^ 9){1 + 6 sin"^ 9) . (31) 

Figure ^ shows the variation of the mass efflux with lati- 
tude. It vanishes at the equator only for /i = — 1. For a par- 
ticular 6, if < 3/ {26 + 1) the maximum of m occurs for 
9 = 90°, while if l/zl > S/{2S+1) it occurs for < < 90°. 
This is shown in Fig. |^(b) for 6 = 4. We simply note in 
passing that the first of these angular dependences of the 
mass loss is reminiscent of some observed intense mass loss 
rates that are thought to occur through equatorial stellar 
winds {e.g. from Be stars observed equator-on). Also, mea- 
surements of the Lyman alpha emission by the satellites 
Mariner 10, Prognoz, Ulysses and SoHO (Bertaux et al. 
1997) have implied that there is higher mass efflux at the 
equator than at the pole. 

We may next introduce the angular momentum efflux 
l{9) = pVrr^L{9) which is the angular momentum loss 
rate per infinitesimal solid angle at the angle 9, 



i{9) = 



{MlY 



sin'+i ( 



1 + p six?" 9. 



(32) 



The variation of i{9) with latitude is shown in Fig. ^. As 
with the mass efflux, the angular momentum efflux va- 
nishes at the equator for /i = —1. In TT91 the angular 
momentum efflux was assumed a priori, from which the 
simplest possible forms for the azimuthal hydromagnetic 
field were obtained. 

2.5. Hydromagnetic field and density 

The azimuthal components can now be obtained 
Rsm'9 fY^-Y 



V4R,9)=XVo 



y/l + 5sm^'9 Vl 



Ml 



(33) 



e=0.5 e=l e=2 




Fig. 3. Co-latitudinal dependence of the mass efflux ih{9) 
for S = 4: m (a) for /i = — 1, in (b) for e = 1 



B^{R, 9) — XBf) 



sin' 9 fR^/Rl-l 



R V 



(34) 



Note that is maximum at the equator. We recall that 
the radial hydromagnetic field and density are of the form 



Vr{R,9)=VoY{R)^ 



Bo 



1 + fi sin^' 9 



Br{R,9) = + fism^'9, 



p{R, 



Po 
YR^ 



{l + 6sm'^' 9) 



(35) 

(36) 
(37) 



2.6. Balance of forces 

At this stage, we still have to deduce an equation for 
Y{R), together with the form for the pressure p(i?, 9) and 
hence the temperature T{R,9). Going back to the r- and 
^-components of the momentum equation - Eqs.(U) and 
(H) - and assuming that the variables R and 9 separate in 
these equations, we must have 



Q{R, 9) = Qo{R) + Qi{R) sin^' 9, 



(38) 



where Q{R,9) is the dimensionless pressure defined by 
PoV„^ 



p{R,t 



-Q{R,t 



(39) 
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Fig. 4. Co-latitudinal dependence of the angular momen- 
tum efflux l{d): in (a) for /i = — f , in (b) for e = f 



Qo represents tlie spherically symmetric part of the pres- 
sure, while Qi includes the effects of the anisotropy. 
Substitution of Eqs.(|8|), (H) into Eqs.(|) and (|) yields 
the following three equations for Qo{R), Qi{R) and Y{R) 



Qi{R) = - 



eY 
A2 



Y -Y, 
l-Ml 



dQi 
di? 



dQo 
di? 



5v^ 

YR* " 

A2 



2fi dY 
R^dR 



1 - R^/Rt 
l-Ml 



2A2 

Vr 



Y -Y^ 



Ml 



d (\-R^/Rl 



{MlfR'^ dR 

2 dY 
R^dR' 



Ml 



YR^ 



(42) 



in which v is the ratio of the escape speed to the radial 
speed at the base of the outflow 



Vo 



(43) 



To understand the interplay between different forces 
involved in the mechanism of this type of wind, let us 
describe one by one the various terms in the above equa- 
tions. All of them are written so that the pressure gradi- 
ent term is isolated on the left-hand side. Equation ( ^o|) 
represents the force equilibrium across the field lines. On 
the right-hand side of this, the various terms represent, 
respectively, the magnetic pressure, the centrifugal force 



and the magnetic tension. The last two equations ex- 
press the equilibrium of forces along the radial direction. 
Equation ( p| ) shows the anisotropic terms, namely, the 
anisotropic part of the gravitational and incrtial forces, the 
centrifugal force and the magnetic tension term, respec- 
tively. The isotropic terms are shown in Eq. (p2|) . These 
are related to the gravitational and inertial forces, respec- 
tively. Equations (^ and (^ can be combined to give a 
single expression for Y{R) 



dr 

dR 
where 

F{R) 



FjR) 
G(i?)' 



(44) 



Y 



YR'^ {MlYR^ e RMl{l-MlY 



(1 + e)Ml - e R'^ 



G{R) 



Ml 
MMlfY 



Rt 



((2- 



e)Ml 



(1 + e)) 



(45) 



2Ml 



Ml 

- 



<1 



Ml)2 



Ml 



Ri 



1 



(46) 



The above differential equation for Y{R) requires a bound- 
ary condition. For convenience, let us choose Y{1) — 1, 
which defines = Vr{R = 1,9 = Q) from Eq.(|5|). /From 
the classical ideal gas law we can now express temperature 
as 



T(i?,( 



mViYR^ 



2fcs(l-h(5sin^ 



-Q{R,( 



(47) 



3. Parametric study of the solution 

(40) 3.1. Critical points 

In order to determine the full solution we need to solve 
Eq.(||), (H), dl^) for Y{R). This is a first-order non- 
linear differential equation which can be integrated nume- 
rically, using a standard routine. A first inspection shows 

(41) that this equation has a singular point at i? = i?* , y = , 
where Ma = 1, or in other words where the radial velo- 
city equals the radial Alfvenic velocity Va = Br/ y/^'np. 
The term singular point is used here in the sense that 
both the numerator and denominator must vanish there. 
In a more general geometry with meridional components, 
this point, known as the Alfvenic point, corresponds to 
the location where the poloidal speed equals the poloidal 
Alfven speed. It is present in magnetic wind models (c./. 
Weber & Davis 1967, Mestel 1968) and is a consequence 
of the steady-state assumption. It delimits the magneti- 
cally dominated region beyond which the torque exerted 
by the magnetic field ceases to dominate over the angular 
momentum carried by the fluid. The exact position of this 
point cannot be found analytically in this case, which com- 
plicates any type of numerical treatment as we shall see 
later. In a similar problem, TT91 find that expanding up 
to fourth order around (i?*, F,), all slopes are allowed. We 
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are thus in the presence of an improper node or star-type 
point. This means that no particular solution is filtered 
out by the presence of this singular point. 



Table 1. Location of both singular points for A = 0.5, 
u = 120, (5 = 4: in (a) for M\ = 0.1 and in (b) for 
Ml = 0.01 




Fig. 5. Topology of the radial dependence of the radial 
velocity Y{R) for a solar-type highly magnetized star with 
A = 0.5, V = 120, Ml ^ 0.01, (5 = 4, e = 1 and ^ = -0.1. 
Note the presence of both critical points 



There is a second singular point, say at {Rx,Yx), 
found by satisfying simultaneously F = G = in Eq. (Q) . 
A first-order analysis around this point shows that only 
two slopes are allowed, one positive and one negative, giv- 
ing rise to a saddle or AT— type point. Thus, this point is 
also a critical point in the sense that it chooses a critical 
solution. In a recent work, Tsinganos et al. (1996) have 
shown that, in general, at the critical points of axisym- 
metric and self-similar solutions, the component of the 
flow velocity perpendicular to the directions of axisym- 
metry (the (/(-direction in this model) and self-similarity 
(the ^-direction in this model), equals the characteristic 
slow/fast MHD wave speed in that direction. 



C^V^ — 



(48) 



with Cs the sound speed and Va , Va,r the total and radial 
components of the Alfven speeds, respectively. The above 
equation is the equation for the speed of an MHD wave 
propagating in the r-dircction. 

Thus, at the critical point where Eq. (^) is satisfied, 
the r-component of the flow speed equals to the slow or 
the fast MHD mode wave speed in that direction. Now, 
in low magnetisation plasmas and in the case of paral- 
lel wave propagation, the slow mode wave speed coincides 
with the Alfven speed. On the other hand, in high magne- 
tization plasmas and in the case of parallel wave propaga- 
tion, it is the fast mode wave speed that coincides with the 
Alfven speed. In the present case we have a degeneracy in 
the sense that there are only two critical surfaces in view 
of the assumption that the critical surfaces are spherical 
(which seems to be a good assumption, at least in the case 
of the solar wind - Exarhos & Moussas 2000). Note that 
due to the non-polytropic assumption, the sound speed 
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(b) 



is ill-defined and we cannot easily determine whether we 
are in a low- or high-magnetisation regime. We conjecture 
that here we have a situation analoguous to low magneti- 
sation plasmas and at the X-type critical point the super- 
Alfvenic radial flow speed equals to the radial component 
of the fast MHD mode wave speed, while the slow critical 
transition coincides everywhere with the Alfvenic transi- 
tion because of self-similarity. 

We should stress at this stage that both the Alfvenic 
point (i?*,y*) and this X— type point are (loosely) called 
equilibrium or critical points since they satisfy the regular- 
ity condition F{R) = G{R) = 0. However, since the first 
one does not filter any solution (i.e. all slopes are allowed 
through it), while the second one (as for any saddle point) 
selects a particular solution sometimes referred to as the 
critical solution, only the X-type singularity corresponds 
to a true critical point. 

We performed a numerical integration of Eqs.(44|), 
15|), (^). The general topology is shown in Fig. | for 
a highly magnetized medium and for representative val- 
ues of fi and e. The position of the critical points was 
found by an iterative procedure. The integration yields a 
single solution crossing the X— type point with positive 
slope and satisfying the boundary condition 1^(1) — 1. 
For this critical solution, the flow starts near the star 
with low speeds and connects to large distances where 
it attains large super-Alfvcnic velocities. There is another 
critical solution, which crosses the X-type point with neg- 
ative slope and is always decelerating. There is a limiting 
value of e above which we couldn't find a wind solution 
satisfying the boundary condition Y{1) = 1. For example, 
if A = 0.5, = 120 and (5 = 4, there is no solution for 
e > 5, if A/^ = 0.1 and for e > 2, if A/^ = 0.01. 

Once the positions of both i?, and Rx are known, 
all other types of solutions can be easily found. They in- 
clude breeze-type solutions that cross the Alfvenic point 
but not the X-type point and reach sub- Alfvenic asymp- 
totic speeds for large distances. 

Table 1 shows the positions of both singular points for 
= 0.1 and = 0.01 and different values of e and /i. 
We have taken A = 0.5, u = 120 and (5 = 4. 
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3.2. Radial velocity 

The initial acceleration decreases as decreases, Figs, 
^(a, b), and also as A increases. In other words, the more 
the outflow is magnetised and rotating, the smaller is the 
initial acceleration and both the magnetic field and rota- 
tion are inhibiters of the initial acceleration of the outflow. 
This may be understood using a simple argument, as fol- 
lows. Assume without loss of generality that 5 = 0, such 



(a) 



H=-0.01 



1000 r 



100 



10000 



£ 100 ^ 




100 



Fig. 6. Radial dependence of the dimensionless radial ve- 
locity Y{R) for a solar-type star with A = 0.5, = 120, 



(5 = 4, e = 1: in (a) for 



0.01, in (b) for, 



0.1. 



that the weight of a parcel of plasma is the same at any 
angle 6, at some fixed radial distance R. Also for sim- 
plicity assume that ii — —1 and e = 1. Then, because of 
the sinusoidal dependence of the radial flow speed with 
the latitude, at the equator we have static conditions such 
that the weight of the plasma is balanced there mainly 
by the centrifugal force and the radial component of the 
Lorenz force (the centrifugal force is always positive, i.e., 
outwards, while in most of the cases examined the radial 
component of the magnetic force is also positive) . Then, 
moving at the pole and at the same radial distance R, we 
find the same inwards gravitational plasma weight (since 
(5 = 0), but no centrifugal and radial magnetic forces act- 
ing there due to their sin^O dependence. The net result is 
a decelerating force, relatively to the cases = A = 0. 



For a mildly magnetized case (Af^ = 0.1), Fig. ^(b) shows 
that the initial acceleration is very high and the radial 
velocity rapidly attains its asymptotic value. 

The asymptotic form of Y{R) for large R can be ob- 
tained after integrating Eq.(B3) 



1 

1 + - 
e 



A2 



InR. 



(49) 



For ^ 1, the asymptotic speed increases logarithmi- 
cally with R, and Y increases with a decrease of or 
a decrease of e and When > 1, the logarithmic 
derivative dF^/dlni?, for large R, is very small (behaving 
like (M^)~^), and Y can be taken as constant. The flow 
is hydrodynamically dominated from R — 1, while both 
singular points, i?* and Rx, almost coincide in the region 
R < 1 {i.e., below the base of the wind). We recall that in 
the hydrodynamic case (Paper I) , Y approaches rapidly a 
constant value, independent of e. 

The cautious reader may have already noted that in 
Eq. (|4^) the asymptotic radial speed is the product of 
a logarithmically diverging isothermal wind speed which 
increases like (lni?)i/3 (see TT91) and Michel's charac- 
teristic speed (rPFg/Aiy/^, where Fb and AI are the 
magnetic and mass fluxes (MacGregor 1996). This is also 
checked by the fact that the effective polytropic index 
7(i? — > oo) — > 1, as seen in Fig. ^. 

3.3. Temperature and effective polytropic index 

In this study by not using a polytropic relationship be- 
tween pressure and density (p oc p^, with 7 the constant 
polytropic index) we avoided constraining the exchange of 
energy so as to keep 7 constant. In our approach, we can 
define an effective polytropic index given by 



7 



dlnp 
dlnp 



A=const. 



(50) 



for each field fine A{e) = const. (Weber 1970, TT91). 
This effective polytropic index is no longer a constant, 
but instead a function of R. The form of 7(i?) is closely 
associated with the variation of temperature with distance 
from the central object. In particular, if 7 < 1 there is 
intense heating and the temperature increases, whereas if 
7 > 1 there is a depletion of heating and the temperature 
decreases (see Fig. 

In Fig. ^ we have chosen to illustrate the two- 
dimensional character of the solution deduced in this pa- 
per, with parameters corresponding to a mildly magne- 
tized object {M\ = 0.1), by plotting the latitudinal de- 
pendence of the radial velocity and density. The velocity, 
is higher at the polar axis while the density is maximum 
at the equator. The angular dependence of these quanti- 
ties is similar to the latitudinal distribution of the density 
and velocity in coronal holes. 
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Fig. 7. Radial clcpcndence of: (a) the dimensionless tem- 
perature at the pole Tp{R)/To, (b) the polytropic index 
7(i?), both for A = 0.5, u = 120, M'i = 0.1, S = 4, e = 1 
and two values of /j.. 



4. Discussion 

We have explored in this paper a new set of exact solutions 
of the steady axisymmetric MHD equations relevant to 
stellar wind problems. 

Analogously to the case with no magnetic field (Paper 
I), wc were able to obtain the latitudinal dependences 
of the different variables and these involve only the 
anisotropy parameters, 5, e and /z. The parameter S is 
associated with the ratio of equatorial to polar density, 
e controls the width of the density and velocity profiles 
for a fixed variation between pole and equator, while n is 
connected to the latitudinal anisotropy of the distribution 
of radial kinetic energy density. The other three parame- 
ters of the problem, which come about in the solution for 
the radial fimctions, arc A, which controls the strength of 
rotation, i/, a measure of the escape speed from the cen- 
tral object and M° which is related the strength of the 
magnetic pressure. 

The degree of coUimation, together with the anisotropy 
in the density distribution, is effectively controlled hy S, e 
and ji. These solutions could be pertinent to the plasma 
dynamics in polar coronal holes, stellar and extragalac- 
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Fig. 8. Contour plots of: (a) the radial velocity 
VriR,9)/Vo and (b) the density p{R,0)/po, for A = 0.5, 
u = 120, = 0.1, 5 = 4, e = 2, ^ = -1 



tic jets and even star-forming regions. The solution with 
/X = — 1 has zero radial velocity at the equator and a con- 
tinuous radial magnetic field across the equator. A dis- 
continuity and associated current sheet arises for /i 7^ — 1. 
Also, the angular velocity of the roots of the; field lines at 
the surface of the central object is such that it resembles 
qualitatively solar observations only for e > 1. 

The topology of the solutions is controlled by two dis- 
tinct singular points. The first one is the usual Alfvenic 
point and corresponds to the location where the radial ve- 
locity of the outflow equals the radial Alfvenic velocity. 

All solutions pass through this high-order singularity and 
it corresponds to an improper node or star-type point. 

There is a critical point present downstream from this 
Alfvcn transition. It is an AT-type critical (or saddle) 
point which is responsible for filtering a single solution 
corresponding to a vanishing pressure at infinity — the 
wind-type solution. In a similar problem, Tsinganos et al. 
(1996) have shown that at this new critical point, the r- 
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component of the flow speed equals the slow/fast MHD 
mode wave speeds in that direction. 

Because the field/streamlines are constrained to keep 
a helicoidal geometry, the pressure gradient for large dis- 
tances has to be sufficiently large to balance the dominant 
tension force. Thus, in some cases the solution is not valid 
outside a cone around the polar axis. This limitation dis- 
appears for > 0.05. 

There is a drastic change in the nature of the solu- 
tion from high to lower values of M\. The solutions for 
high M\ manifest all the characteristics of a typical hy- 
drodynamic wind, namely very large acceleration at the 
base of the atmosphere and temperature decaying very 
rapidly with distance. The magnetically dominated cases 
(M^ <C 1) show lower acceleration at the base of the at- 
mosphere and an almost isothermal atmosphere at larger 
distances. 

The density anisotropy [5] greatly favours the accel- 
eration of the wind close to the base, while the strength 
of rotation (A) and the gravitational field [v) slows down 
the initial speed. The solution far away from the source 
seems to be quite insensitive to e. The influence of /i, on 
the other hand, is such that a decrease of not only di- 
minishes the relative importance of the magnetic effects 
by decreasing the size of the magnetic lever arm, i?*, but 
also increases the value of the initial acceleration as well 
as the asymptotic radial velocity at large distances. The 
asymptotic plasma temperature increases for lower values 

This simple but self-consistent two-dimensional solu- 
tion of the MHD equations has nevertheless some limi- 
tations. The assumption of separation of variables could 
be dropped, although we should not forget the nonlinear 
character of the governing coupled partial differential 
equations. The assumption of the neglect of the merid- 
ional components of the magnetic field and flow {Vg,Bg) 
which limits their possible applications, could also be elim- 
inated within the framework of the method of separation 
of the variables. Finally, as we've already discussed, the 
present investigation does not take into account the de- 
tailed energetics in the outflow. Nevertheless, the present 
study indicates that inherently non-spherically symmetric 
solutions of the MHD equations which do not have the 
previous limitations may exist. This ambitious undertak- 
ing to flnd them remains a challenge for the future. 
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